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ABSTRACT 

We use Spitzer MIPS data from the FIDEL Legacy Project in the Extended Groth Strip to analyze 
the stellar mass assembly of massive (A / t>10 11 M.q) galaxies at z<2 as a function of structural 
parameters. We find 24 /an emission for more than 85% of the massive galaxies morphologically 
classified as disks, and for more than 57% of the massive systems morphologically classified as spheroids 
at any redshift, with about 8% of sources harboring a bright X-ray and/or infrared emitting AGN. 
More noticeably, ^60% of all compact massive galaxies at z=l-2 are detected at 24 /xm, even when 
rest-frame optical colors reveal that they are dead and evolving passively. For spheroid-like galaxies 
at a given stellar mass, the sizes of MIPS non-detections are smaller by a factor of ^1.2 in comparison 
with IR-bright sources. We find that disk-like massive galaxies present specific SFRs ranging from 
0.04 to 0.2 Gyr" 1 at z<l (SFRs ranging from 1 to 10 M® yr" 1 ), typically a factor of 3-6 higher than 
massive spheroid-like objects in the same redshift range. At z>l, and more pronouncedly at z>1.3, 
the median specific SFRs of the disks and spheroids detected by MIPS are very similar, ranging from 
0.1 to 1 Gyr~ x (SFR=10-200 7W© yr _1 ). We estimate that massive spheroid-like galaxies may have 
doubled (at the most) their stellar mass from star-forming events at z<2: less than 20% mass increase 
at 1.7<z<2.0, up to 40% more at l.l<z<1.7, and less than 20% additional increase at z<l. Disk-like 
galaxies may have tripled (at the most) their stellar mass at z<2 from star formation alone: up to 
~40% mass increase at 1.7<z<2.0, and less than 180% additional increase below z=1.7 occurred at a 
steady rate. 

Subject headings: galaxies: evolution — galaxies: starburst — galaxies: elliptical — galaxies: forma- 
tion — galaxies: photometry — galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

The formation and evolution of massive 
{M>1Q 11 Mq) galaxies is one of the most studied 
topics in extragalactic astronomy during the last decade. 
Early expectations from hierarchical galaxy formation 
models considered that star formation began in low 
mass systems which built more massive galaxies throug h 
sequential merging (iB augh et al.lll996l : ICole et al.ll2000f ). 
in a similar process to the gro wth of structures i n 
Cold Dark Matter simulations ([Springel et al.l l2005f ). 
In apparent contradiction with hierarchical assembly, 
the finding of a substantial population of massive 
galaxies at z>l (lElston et al.l[l9 88; Hu ghes et all 119981 : 
iFranx et aT1l2"00 3: Glazcbrook et al. 2004), some of them 
containing old stellar populations and evolving passively 
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(accor ding to their red optical colors - 4Daddi et al.l 
Redd v et al.l 120051— , and spectra - iKriek et al.l 
2008; Cimatti et al. 20 08)-), seems to fav o r a downsizing 
form a tion scenario (ICowie et al.l 119961: I Heavens et all 
20041: Uuneau et al.1 120051: IPerez-Gonzalez et al.l 120051 : 
Bundv et al.l [2006). This population of z>l massive 
galaxies accounts for a significant fraction of the lo- 
cal stellar mass density (^20% as early as z~2, and 
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~10% at z~4: iFontana et alj 120061 : 1 Arnouts et~alll2007t 
IPerez-Gonzalez et al.l [20081. The discovery of such a 
population reinforced the idea that both stars and their 
host galaxies are coeval (resembling a monolithic-like 
collapse), and consequently, no expectations of structure 
evolution in these galaxies should be expected. For this 
reason, the recent observational evidence showing that 
the most massive gal axies were much more co mpact 
in the past (iDaddi et all 120051: iTruiillo et al.1 120061: 
lLonghetti et all 120071: ICimatti et all l2008f) has been 
surprising, and has again opened the question of how 
the stellar populations of these galaxies were assembled 
into their present shape. 

The size evolution of the most massive objec t s sinc e 
z~2 has been characterized by Truii llo et al.l (|2007f) . 
These authors found that, at a given stellar mass, disk- 
like objects at z~1.5 were a factor of two smaller than 
their present-day counterparts. For spheroid-like ob- 
jects, the evolution has been even stronger: they were 
a factor of four smaller at z~1.5 than nearby similar 
mass ellipticals. In addition, the stellar mass densi- 
ties of these high-z galaxies were almost two orders of 
magnitude higher than objects of the same mass to- 
day. These superdense gal axies have been fo u nd at even 
higher (z~2.5) redshifts flZjr m_et_aL| 120071 : iToft et all 
l2007t Ivan Dokkum et aT1l2008h , adding more controversy 
to the debate about the formation and evolution of mas- 
sive galaxies. 

Two processes have been proposed to allow the su- 
perdense high-z galaxies migrate to the local stellar 
mass-size relation. The first process is dissipation- 
less (with absence of star formation) merging. Given 



the high metal abundances and old ages of the stel- 
lar p opulation present in l ocal massive elliptical galaxies 
(e.g. jGallazzi et aL|[200llSanchez-Blazquez et al.ll2006t 
iJimenez et al.l 20071). these m ergers should preferentially 
be dry ( De Lucia et al.l[2006f ). and occur between z~1.5 
and z=0, the epo ch when the red sequence appears 
(|Labbe et al.ll2007T ). In this context, a particular effective 
size evol utionary mechanism ( t e^^M ) has been pro- 
vided by iBovlan-Kolchin et al.l (2006) through head-on 
mergers of galaxies. The second possibility is the smooth 
envelope accretion scenario (jNaab et al.ll2007i ) , where ac- 
creted stars (mainly provided by minor mergers) form 
an envelope whose size increases smoothly at decreasing 
redshift. 

The goal of this paper is to explore the evolutionary 
paths followed by the most massive galaxies and their 
dependence on the morphology. To do this, we quantify 
the growth in stellar mass via star formation events of 
massive (A^>10 11 M.®) galaxies as a function of size 
and brightness profile shape up to z^2. We base our 
discussion on the characterization of the dust infrared 
(IR) emission of these systems, which is linked to the 
amount of recent star formation and/or the presence of 
obscured AGN. This IR-based study is complementary 
to the more classical approach to the characterization of 
the evolution of massive ellipticals based on rest-frame 
optical properties. 

Throughout this paper, we use a cosmology with Ho = 
70 kms^Mpc" 1 , Q M = 0.3 and A = 0.7. All mag- 
nitudes refer to the AB sy stem. Th e resul ts for stellar 
masses and SFRs assume a Ch abrierl ()2003[ ) initial mass 
function (IMF) with 0.1<oW<100 M®. 

2. DATA DESCRIPTION 

2.1. The sample 

To analyze the star formation properties of the most 
massive galaxies as a function of morphology, we use 
the catalog of 831 A"-band selected massive galax- 
ies (A / J>10 11 M.®) in the Palomar Observ atory Wide- 
Field Infrared (POWIR)/DEEP-2 survey rfDavis et all 
20031: iBundv et ; alJl2006t IConselice et al.ll2008ft for which 
Truiillo et al.1 (|2007h provide redshifts, stellar masses, 
and structural parameters (sizes and Sersic indices). 
These data, jointly with the Spitzer/MIPS fluxes mea- 
sured in the observations carried our by the FIDEL 
Legacy Program in the Extended Groth Strip (EGS), 
allow a detailed analysis of the star formation properties 
of the most massive galaxies as a function of morphology 
up to z^2. 

The sam ple is describe d in d etail in IConselice et al.l 
(J2007J ) and iTruiillo et all (|2007h . Briefly, the K-band 
survey covers 2165 arcmin 2 in the EGS and has a depth 
K=22.5 mag (5cr). Only 710 arcmin 2 are covered simul- 
taneously with HST/ACS v- and i-band imaging from 
the All- Wavelength Ex tended Groth Strip International 
Survey (AEGIS, iDavis et aT] 120071) . so reliable struc- 
tural parameters can only be measured for 831 galaxies 
within the entire POWIR/DEEP-2 survey in EGS. Half 
of those 831 galaxies has spectroscopic redshifts based 
on optical d ata obtained by th e DEEP-2 Galaxy Red - 
shift survey ([Davis et al.ll2003f ). IConselice et~aT1 (|2007l ) 
estimate photometric redshifts for the rest of sources 
with an accuracy Az/(1 + z)=0.025. The 831 galax- 
ies with M>10 n M® in the EGS lie in the red- 



shi ft range 0.2<z<2. Stel lar masses were estimated 
bv lBundv et al.l (I2005L I2006T) using the e x ponen tial star 
formation mod els of iBruzual fc Chariot] (|2003[ ) with a 
I Chabrierl (2003) IMF and various ages, metallicities and 
dust c ontents included. As shown by IConselice et al.l 
(2007), typical uncertainties in the stellar masses are a 
factor of ~2 (typical of any stellar population study; 
see, e.g., |Pcrcz-Gonza lez et al.l [20031 iKa uffman n et al.l 
20031 (Papovich et all I2006L I Font ana et all I2006L and 
Perez- Gonzale z et al.l J20 08V As discussed in detail in 



Conselice et al.l (|2007f ). this factor includes the effects 



of the photometric redshift uncertainties, the errors in- 
herent to solution degeneracies, and the choices of the 
IMF and the ste llar emission library. For example, using 
iMarastonl (2005) models (with an improved treatment of 
the TP-AGB stellar evolution phase) would produce a 
20% (at most) systemat i c dec rease in the mass estima- 
tions. Using a ISalpete"r1 (|1955f ) IMF would increase the 
st ellar masses by a co nstant factor of 0.25dex. 

iTrujillo et al.l (120071 ) estim ated (circularized) half-light 
radius (r e ) and Sersic (1968) indices (n) for all the galax- 
ies in our sample. They used z-band HST/ACS images 
to fit surface brightness profiles and divided the sample 
in disk-like and spheroid-li ke galaxies accord i ng to the 
valu e of the Ser sic index. R avindranath et all (|2004L see 
also Andrcdak is et al.1 11995) demonstrated that nearby 
galaxies with n<2.5 are mostly disks, while spheroids 
are characterized by high Sersic indices, n>2.5. They 
also performed simulations to check that the Sersic index 
obtained from HS T data can be used as a morphology 
indicator at z>0. ITruiillo et al.l (|2007l ) extended these 
simulations to prove that the structural parameters are 
robust against the effects introduced by if -corrections, 
AGN contamination, and surface brightness dimming. 

For our sample, visual inspection of the ACS i-band 
images by one of the co-authors (I.T.) was used to clas- 
sify the sample in 4 types: ellipticals/lenticulars, spirals, 
irregulars, and mergers. Comparing this visual classifi- 
cation with the one based on Sersic indices, we find that 
the visually confirmed spheroids present <n>=4.8±1.5, 
and the rest of sources have <n>=1.7±1.7. There is 
a 6% contamination of visually identified spheroids in 
the n<2.5 sample, and a 20% contamination of visu- 
ally identified disks in the n>2.5 sa mple (comparable to 
the 5 % and 19% contaminations in iRavindranath et al.l 
2004). The fraction of interloper disks decreases to 7% 
at n>4 and 4% at n>5. A 4% contamination is typi- 
cal of other works based on visu al or quantitative mor - 
ph ological class if icatio ns such as IConselice et all (|2007l ) 
or IBundv et al.l (|2005| ). The visual test shows that the 
galaxies with n>4 form a robust (almost uncontami- 
nated) sample of spheroid-like sources, and n<2.5 galax- 
ies are mostly disks. Galaxies with 2.5<n<4 are most 
probably spheroid-like galaxies with some contamination 
of SO galaxies and early-type spirals. 

2.2. UV-to-MIR photometric properties of the sample 

The positions of the 831 massive galaxies in the EGS 
were cross-correlated (using a 1" search radius) with our 
own reduction and catalogs of the Spitzer IRAC survey 
of the EGS. Using the same sim ulation method described 
in lPerez-Gonzalez et al.1 (|2008t ). we found that this cata- 
log is 75% complete at ~1.5 ^Jy ([3.6]=23.5 mag), which 
corresponds to 8<r detections. Our IRAC photometry is 
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Fig. 1. — Postage stamps and SEDs for four typical MlPS-detected galaxies in our sample of massive galaxies. Left panels show 
10"xl0" RGB composite images built from HST/ACS V and i frames. In the middle column, MIPS 24 (im images of size 40"x40" are 
depicted, with the red square showing the area covered by the ACS postage stamp. In all images, North is up and East is left. The right 
columns show the SEDs of each galaxy, fitted to stellar population and dust emission models which are used to estimate photometric 
redshifts, stellar masses and SFRs (these parameters are given in each SED plot). The two upper rows show examples of disk-like galaxies: 
EGS142126.97+531137.4, a galaxy at z=0.67; and EGS142013. 18+525925.0, lying at z=1.65. The two lower rows show examples of 
spheroid-like galaxies: EGS142021.47+525543.4, a galaxy at z=0.63, and EGS142125.76+531622.8, placed at z=1.70. 



consistent with that published by iBarmbv et al.l (|20Q8T ) 
for the same dataset (but their own reduction and cat- 
aloging) within 0.1 mag (typical absolute uncertainty of 
IRAC fluxes) for 75% of the sample, and within 1-er error 
for virtually all sources. 

We found IRAC counterparts down to [3.6]=23 mag 
for all of the 831 galaxies in Trujillo's sample. For 
151 sources (18% of the total), the IRAC sources were 
blended with nearby objects, but still resolved (the 
separation was larger than 1"). As described in de- 
tail in Appendix A of IPerez-Gonzalez et al.l (2008), for 
these sources we obtained multi-wavelength photome- 
try (including Spitzer fluxes) using a dcblending al- 
gorithm based on the deconvolution of the IRAC and 
MIPS images. The method takes the known positions of 
the blended sources obtained from optical/NIR ground- 
based images and the PSFs for the different images and 
obtain separated flu xes for the blended sources (see also 
IGrazianet al. 2006). The method relies on the moder- 
ate resolution of the IRAC images (~2" FWHM, not that 
different from an optical ground-based image, but very 
stable) , which allows the deblending of sources separated 
by more than -1" (half of the FWHM). For the MIPS 
24 /im images, the resolution is worse (~6" FWHM) but 
the IRAC data can be used to assign most probable coun- 
terparts and help with the deblending. In any case, the 
main results in this paper (the MIPS detection fraction 
and the statistics of specific SFRs) remain virtually un- 
changed (less than 5% random changes at all redshifts) 
when we remove the 151 sources with blending problems. 

We measured consistent aperture photometry in sev- 
eral UV, opti cal, NIR, and MIR bands wi th the method 
described in IPerez-Gonzalez et al.l (|2008l ). The multi- 
wavele ngth dataset is outlined briefly in iVillar et~aT1 
(2008) and will be characterized in detail in Barro et 
al. (2008, in preparation). More noticeably, our merged 
photometric catalog includes MIPS fluxes at 24 /im ob- 
tained from aperture photometry in the GTO and FIDEL 
survey (DR2) data in the AE GIS/EGS field (jDavis et all 
l2007t ISvmeoni dis et al. 2007). Following the same pro - 
cedure described in IPerez-Gonzalez et al.l (|2005, 2008), 
we used the DAOPHOT software package in IRAF 5 to 
detect sources (using a 3cr detection cut above the lo- 
cal sky noise) in the MIPS images and measure aperture 
photometry with a PSF fitting method. 

Based on simulations consisting on the creation and re- 
covery of artificial sources in these images, we found that 
our 24 fim catal og of the EGS is 75% c om plete at 35 /zJy 
(cons istent with iPapovich et aill2004l and lTreister et all 
2006). i^(24)=35 fiJy corresponds to ~6cr detections 
for the average sky noise in the FIDEL DR2 images. 
Our 24 fim catalogs are cut to 3a detections, which 
translate to the range 14-17 /tJy, depending on the lo- 
cation on the sky due to small differences in exposure 
time and the effect of confusion (both presenting small 
spatial variations througho ut the image). As done in 
IPerez-Gonzalez et al.1 (J2005), we tested the reliability of 
the MIPS detections by analyzing the probability of hav- 
ing a counterpart within the search radius (1") in other 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc., under cooperative agree- 
ment with the National Science Foundation. 



optical/NIR (ground-based and IRAC) bands for a ran- 
dom position on the sky of the 24 fim image. Having a 
counterpart in 3 different bands within the EGS dataset 
has a negligible probability (1.4%), so we conclude that 
(virtually) all the MIPS detections within our sample are 
not spurious. Figure [1] shows postage stamps and spec- 
tral energy distributions (SEDs) of two typical examples 
of disk-like galaxies and two spheroid-like galaxies de- 
tected by MIPS at 24 /im (also one 70 /im detection 
included). 

Using the measured MIPS 24 /im fluxes, the esti- 
mated 280 nm synthetic fluxes inter/extrapolated in the 
spectral energy distribution fits, and the spectroscopic 
and p hotometric redshifts published by iTruiillo et al.l 
(2007), wc obtained total (unobscured plus obscured) 
SFRs for each galax y in t he same way explained i n 
IPerez-Gonzalez et al.l (|2008l see also iBell et all I2005T) . 
Briefly, the MIR fluxes at rest-frame wavelengths longer 
than 5 /im are fitted to dust emission models (from sev- 
eral libraries) and the IR-based SFRs are obt ained from 
integrated total IR (using the cal ibration in | KennicuttJ 
1998) and rest-frame 24 /im (see lAlonso-Herrero et all 
2006) luminosities from the fits (averaged through all 
template libraries). The IR-based (obscured) SFR is 
then added to the UV-based (unobscured) SFR to ob- 
tain the total SFR. As discussed in that paper, total 
SFR estimates should be good within a factor of 2. The 
SFRs discu ssed in the follow ing Sections were estimated 
assuming a IChabrieil (J20 03f) IMF, obtained by dividing 
the re sults obtai ned with the c alibrations in iKennicut/n 
(119981 valid for a lSalpeterlll955l IMF) by a 1.8 factor. 

3. RESULTS AND DISCUSSION 

3.1. IR emission of the most massive galaxies at z<2 

Figure [5] shows the location of the MIPS detections in a 
stellar mass-size diagram for massive galaxies as a func- 
tion of the concentration index. Out of the 831 galax- 
ies in Trujillo's sample of massive galaxies, 485 sources 
(58%) w ere clas s ified a s spheroid-like based on their high 
(ra>2.5) ISersid (|1968l ) indices, and 346 (42%) as disk- 
like sources (n<2.5). Among the disk-like sources, 322 
(93% of all disks) are detected by MIPS at 24 /im with 
a minimum flux i 7 '(24)=15 /iJy, and 137 galaxies are de- 
tected at 70 /im down to F(70)=0.5 mJy. Among the 
spheroids, 297 galaxies (61% of the total) are detected 
at 24 fim down to F(24)=14 /iJy, and 84 sources are de- 
tected at 70 fim down to -F(70)=0.7 mJy. If we consider 
only the MIPS 5a detections (i.e., more statistically reli- 
able sources, although our simulations reveal that all our 
3cr 24 fim sources are reliable to the 99% level; see Sec- 
tion E2]), the MIPS detection fractions decrease to 92% 
for disky systems and 52% for spheroid-like sources. 

Table Q] shows the total number of sources and MIPS 
detection fractions (for fluxes above the 3cr level) as a 
function of redshift and morphology. We consider the 
results based on the morphological classification using 
the Sersic indices and the direct visual inspection of the 
images. The MIPS detection fractions for n<2.5 galax- 
ies and visually identified disks are almost identical. Vi- 
sually confirmed spheroids present less MIPS detections 
than the n>2.5 sources, although the difference is small 
(<10%) and consistent with the 20% contamination of 
visually identified disks in the n>2.5 sample, most of 
them having 2.5<n<4.0 (see Figure [3]). 
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Fig. 2. — Stellar mass— size distribution for different redshift bins of our massive galaxies separated in disk-like (left panels) and spheroid- 
like (right panels) types. Galaxies detected by MIPS at 24 /im are plotted with filled stars, while open stars show MIPS non-detections. 
Red symbols are galaxies whose MIPS emission is identified with obscured star formation, and green symbols depict galaxies which harbor 
an X-ray or/and IR-emitting AGN. 



TABLE 1 

MIPS DETECTION FRACTION AND SPECIFIC SFR STATISTICS AS A FUNCTION OF MORPHOLOGY. 
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Figure shows that there is basically no difference 
between the loci occupied by MIPS detected and un- 
detected galaxies in the stellar mass-size plane. How- 
ever, for spheroid-like objects at a given stellar mass, 
MIPS non-detections are sm aller than IR-brig ht sources 
by a factor of ^1.2 (see also lZirm et al.l 12007 results at 
higher z). This suggests that early (i.e., z>2) massive 
star formation events left even more compact remnants 
than starbursts taking place at z<2, maybe reflecting the 
higher density conditions of the primitive Universe. 

Figure [3] shows the MIPS 24 /mi detection fractions as 
a function of structural parameters. This Figure confirms 
the bias of the n<2.5 sample towards galaxies with on- 
going (possibly extended through the disk) star forma- 
tion or harboring an IR-emitting AGN. Figure [3] demon- 
strates that virtually all (80-90%) the n<2.5 galaxies are 
detected by MIPS at all redshifts and (almost) indepen- 
dently of the size of the galaxy. In contrast, the spheroid 
sample is biased towards more quiescent systems. There 
is still a non-negligible fraction (6%) of galaxies classi- 
fied as disky which fall below the MIPS detection limit 
or do not present any IR emission, i.e, they have low-level 
star formation, no dust, or are completely quiescent. Ac- 
cording to Figure [31 most of them lie at z>l (59% of all 
disk-like non-detections) and tend to have comparatively 



smaller sizes: the MIPS detection fraction decreases from 
90% for disk-like galaxies with r e >4 kpc to 70-80% for 
r e <1.5 kpc systems. All these sources present very red 
SEDs (see Figure gj). 

On the contrary, the spheroid-like sample includes at 
least a 63% of "active" galaxies. Moreover, some more 
MIPS undetected spheroids may have some star forma- 
tion activity or harbor an AGN, since some of the SEDs 
in the upper-left panel of Figure [4] present a signifi- 
cant emission in the UV, probably arising from young 
stars. Most of these UV-bright galaxies lie at z>l, and 
the MIPS 24 /im flux upper limits 6 are consistent with 
the MIR emission from a typical Sc galaxy. Figure [3] 
shows that MIPS detections are more common among the 
largest spheroid-like galaxies, especially at z>l. Interest- 
ingly, at z>l the MIPS detection fraction stays roughly 
constant at r e <3 kpc. It may even increase (up to 70%) 
for very compact (r e <l kpc) spheroid-like galaxies, al- 
though there are 3 caveats to this result: 1) the number 
of sources with r e <l kpc is small (~30), so the uncertain- 
ties in these bins are of the order of 20-30%; 2) for z>l 
and r e <l kpc we are reaching the resolution limit in the 

6 The flux upper limits at 70 /im have been omitted from Fig- 
ure |4] for clarity, given that very few sources are detected at this 
wavelength. 
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HST/ACS images, and consequently the size measure- 
ments count with large uncertainties; and 3) these high- 
z compact galaxies may be dominated by an obscured 
AGN (since the galaxies are detected by MIPS a nd show 
no X-ray emission; see also iTruiillo et al.l I2007I ) , which 
may bias the size measurements. 

Our results ma y be co mpared with those published by 
iRodighiero eTal! (|2007f) . They found that 20% of the 
most securely identified spheroids 7 at 0.3<z<1.0 are de- 
tected during phases of prominent activity. The evidence 
is the detection at 24 fim. above 80 fiiy or at radio wave- 
lengths above 40 /iJy. Our MIPS data are deeper than 
theirs, but if we cut our catalogs to the same flux limit 
(see Figure [3]), we obtain an average ~20% detection 
fraction at z<l for the most concentrated objects with 
n>4, probably well correlated with their sample of bona 
fide spheroids, although our selection only includes the 
most massive galaxies (A4>10 n Mq), and Rodighiero's 
selection is not based on stellar mass. 

If we consider the high-z (z>l) galaxies in our sam- 
ple, our re sults are also consistent with those found in 
IPapovich et~aT] (|2006l ) for Distant Red Galaxies (DRGs). 
These galaxies have a typical stella r mass .M^IO 1 * Mq 
and a mean redshift z~2 (see also iGrazian et all l2007t 
iPerez-Gonzalez eT al. 2008). Given that our sample has 
a stellar mass cut M>10 Mq and a redshift cut at 
z=2, the low-redshift, high-mass tail of the general DRG 
population must be included in our selection. Indeed, 



7 Note that Rodighicro ct al. (2007) c ut their sampl e to cl ear E 
and E/S0 galaxies as classified visually by[Bundy ct al. (2005|), but 
do not include Bundy's SO type in their analysis. 



we have 27 DRGs in our sample, with an a verage red- 
shift <z> = 1.50±0.22. IPapovich et~all (j2006h found that 
roughly 50% of DRGs are detected by MIPS at 24 /j,m 
down to 80 /zJy. We find a ^60% detection fraction for 
massive galaxies in the low redshift tail of DRGs (l<z<2) 
with deeper data (75% completeness level at 35 /xJy). 
Among the 27 DRGs in our sample, 25 (93%) are de- 
tected by MIPS. This detection fraction is higher than 
the average f or l<z<2 galax i es, bu t still consistent with 
the results in IPapovich et al.1 (|2006f ). who argue that z<2 
DRGs are mostly heavily extincted starbursts part of the 
class of dusty EROs at z> 1, and find a ~75% MIPS de- 
tection fraction for their M>10 n Mq DRGs at z<2. 
Indeed, 18 DRGs in our sample are EROs (16 MIPS de- 
tections) , and we count with a total of 305 EROs in our 
entire sample. 

3.2. Spectral energy distributions 

Figure [4] shows the SEDs of all the massive galaxies in 
our EGS sample divided into morphological and activ- 
ity types. The two panels on the r ight show the SEDs 
for galaxies with an X-ray detection (B armbv et al.ll2006l 
iNandra et alj2007t see also lConselice et al . 2007) or clas- 
sified as IRAC power-law galaxies (PLGs, iDonlev et al.l 
2007i), i.e., sources which most probably harbor an X-ray 
and/or IR emitting AGN. There are 68 X-ray emitters, 
60 of them with MIPS detection, and 5 PLGs in total, all 
of them with X-ray emission and 4 with MIPS emission. 

The SED distribution of the spheroids present a lower 
scatter in the UV/optical, being very similar to a spectral 
template of an elliptical. In contrast, there is a very pop- 
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sources in each type for the spheroid-like and disk-like samples. 



ulated tail of disk-like galaxies with UV/optical fluxes 
brighter than a template for a typical Sc galaxy, most 
probably linked to a recent star burst. The MIR emis- 
sion is consistent with the PAH spectrum of a late-type 
spiral galaxy (see also Figure [T]), but can b e as high as 6- 
10 tim es the flux of the Sc template from iPolletta et al.l 
(2007|. Spheroid-like galaxies with a MIPS detection 
present a lower 24 /mi median flux (82 /iJy, with the 
quartiles being 42 //Jy and 188 /tJy, and the average 
160 /tJy) than MIPS disky sources (190^ M Jy, and the 
average 250 /iJy). Most of the sources identified as AGN 
present relatively bright fluxes at rest-frame wavelengths 
between 2 and 10 /mi, revealing the presence of very 
hot dust heated by the central supermassive black hole 
and emitting in the NIR and MIR. In several cases, this 
NIR/MIR emission hides the 1.6 /mi bump, typically 
seen in galaxies whose spectrum is dominated by stars 
rather than dust. 

3.3. Specific Star Formation Rates 

Figure [5] shows the specific SFRs of massive galaxies 
as a function of redshift and morphology. Table [1] gives 
the median and quartiles for different redshift ranges and 
morphological types (obtained from Sersic index and vi- 
sual classifications). The median specific SFRs increase 
by less than O.ldex and 0.02dex for spheroids and disks, 
respectively, when considering MIPS 24 /mi detected 
above the 5<7 level. These increments do not affect the 



following results significantly. 

When segregating the sample based on the Sersic in- 
dices, we find that the specific SFRs of spheroid galaxies 
evolve as (l-fz) 55±0 - 6 from z=0 to z=2, while the evo- 
lution for disk-like galaxies goes as (l+z) 36±0 - 3 . If we 
consider the results based on the visual classification, the 
evolution is more pronounced for spheroids and almost 
identical for disks: (l+z) 6 - 4±08 evolution for the former, 
and (l+z) 3 - 4±0 - 2 for the latter. 

The specific SFRs used in Figure [5] have been esti- 
mated by adding the unobscured SFRs obtained from 
UV data (at rest-frame 280 nm) and the obscured SFRs 
from IR data (using the total IR lu minosity) as ex- 
plained in iPerez-Gonzalez et al.l (|2008l ). The ratio be- 
tween these two quantities allow the estimation of the 
global obscuration of the recent star formation in each 
galaxy (only for those detected by MIPS). On average, we 
find that extinctions for MIPS detected galaxies increase 
with redshift, ranging from <^4(V r )> = 1.0±0.5 mag 
at z<0.5 to <^4(F)> = 1.5±0.6 mag at z^l and 
<A(y)>=2.0±0.7 mag at z^2. These values are consis- 
tent with t ypical attenuations foun d for IR-bright galax- 
ies by, e.g.. lRigopoulou et al.1 (|2000h . and the evolution in 
th e extinction propert ies of the UV SF R density found 
by iTresse et al.l (|2007fl . According to ISvmeonidis et all 
(2007), even larger extinctions (up to a factor of ~100) 
are needed to match SFRs obtained from IR or radio data 
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Fig. 5. — Specific SFRs as a function of rcdshift and morphology 
(for galaxies not identified as bright AGN). Galaxies detected at 
24 fim are plotted with open (disks) and filled (spheroids) black cir- 
cles, while gray symbols show upper limits for sources not detected 
by MIPS. Red and blue crosses represent the median and quartiles 
for the distribution of specific SFRs in the different redshift ranges 
used in Truiillo ct al. (2007, red widest lines referring to spheroids 
and blue narrowest to disky galaxies). Green curves show the ex- 
pected positions of galaxies which would multiply their stellar mass 
by 5/4, 2, and 4 between their redshift and z=0 if they maintained 
a constant SFR. Horizontal dashed lines show constant SFR values 
for the median stellar mass of our sample (1.6X10 11 Mq)- 



and SFRs obtained from [Oil] spectroscopic observations 
of IR-bright sources in the EGS. 

Figure \5\ shows that below z=l.l spheroid-like galax- 
ies present very low specific SFRs. On average, they 
would increase their stellar mass by less than 25% at 
0<z<l if they maintained a constant SFR. The global 
mass increase (in the form of newly-formed stars) for all 
spheroid-like galaxies is less than 10% if we take into ac- 
count the 38% of z<l spheroids which are not detected 
by MIPS, and even lower (—5%) if we only consider the 
visually identified spheroids. In contrast, disk-like galax- 
ies could typically double their mass from z=l to z=0 
due to newly formed stars if they maintained a constant 
SFR, with little change due to the very few galaxies (less 
than 5%) for which we only have SFR upper limits. 

In practice, IR-bright i ntense star-forming burs ts are 
not expected to last long (|Mihos fc Hernq uist 1994), so a 
galaxy most probably will not maintain a high SFR level 
for several Gyr. The higher SFR values are expected to 
be maintained for shorter periods, since gas exhaustion 
and supernova winds (and e ven AGN activity) w i ll help 
to suppress star formation. iMihos fc Hernquist] (1996) 
simulations of mergers predict the triggering of a very 
intense and short starburst event (probably detectable in 
the MIR by MIPS) lasting a few tens of Myr (and occur- 
ring in late stages of the merger, when the galaxies are ac- 
tually joining) for encounters of galaxies with an already 
formed bulge. Encounters of disky galaxies would trig- 
ger less intense bursts lasting longer (100-200 Myr) and 
occurring earlier in the merger process. For the observed 
specific SFRs in our sample, those short and intense star- 



bursts would add up less than 0.01% (for each merging 
event) to the total stellar mass of a typical spheroid-like 
galaxy at z<l. This very small fraction of young stars 
would be hidden by the predominant old stellar popu- 
lation and be undetectable in local ellipticals. For the 
disky galaxies at z<l, which present specific SFRs as 
high as 0.2 Gyr -1 , the burst strength (ratio of the newly 
formed stars to the global stellar mass) could be as high 
as a few percent (for each merging e vent), typical for 
star-f orming galaxi es at low-reds hifts (|Kauffmann et all 
I200& iPerez-Gonzalez et al.ll2003l) . 

At z>l, the specific SFRs of massive galaxies are 
higher than 0.1 Gyr -1 , both the active spheroid-like 
(note that there are 40% of spheroids which are not de- 
tected by MIPS) and disky systems are forming stars at 
approximately the same rate, and the number of quies- 
cent galaxies (those not detected by MIPS) is less than 
—50% for both types. It is interesting to notice that most 
galaxies (disks and spheroids) have significant amounts 
of dust, since they are detected at 24 /im. If some of 
these galaxies are the progenitors of nearby ellipticals, 
that dust should have disappeared somehow or it is now 
very cold and may only be detected at very long wave- 
lengths (A> 100-200 /im) and low fluxes. 

For typical burst durations, and even for star- forming 
events with a constant SFR and lasting up to 1 Gyr, 
the maximum increase in stellar mass would be —15% 
at l.Kz<1.4, -25% at 1.4<z<1.7, and -50% at 
1.7<z<2.0, for both spheroids and disks. This means 
that a significant fraction (more than 50%) of the stel- 
lar mass of z>l massive gal axies was assembled at z>2 
dPerez-Gonzalez et alj|2008t see also iBauer et al.ll2005l 
iFeulner et al.ll2005l iPapovich et alJl2006f ). Moreover, we 
find that about 40% of spheroids at z— 1.8 are almost 
"dead" (they present low SFR levels based on IR and 
UV data) and evolving passively, or may be experiencing 
a quiescent period. Note that most spheroid-like galaxies 
at z>l would be qualified as passive based on optical col- 
ors alone, but the MIPS data reveals that —50% of them 
are experiencing dusty starbursts and 10% more harbor 
(also) obscured AGN. 

We can estimate how much stellar mass galaxies typ- 
ically assemble through star formation from z=2 to the 
present (i.e., the star formation efficiency to increase the 
mass of a galaxy in the last 10 Gyr) if a galaxy follows 
the specific SFR evolution depicted in Figure [5] We as- 
sume that the SFRs remain constant within each redshift 
interval; since the starbursts probably last 50-200 Myr, 
as discussed earlier, the following figures would be an 
upper limit. Adding all the mass formed from z=2 to 
z— 0, we estimate that a disk-like galaxy could increase 
its stellar mass by up to a factor of 3.2±0.5 in the last 
10 Gyr: 1.4 times increase at 1.7<z<2.0 and an almost 
constant 10-20% increase in each of our 5 redshift inter- 
vals at z<1.7. A spheroid-like galaxy could increase its 
stellar mass by up to a factor of 1.8±0.3: 1.2 times in- 
crease at 1.7<z<2.0, 10-20% in each of our two intervals 
at l.l<z<1.7 and less than 5% in each of the three inter- 
vals at z<l.l. These figures are almost unchanged (<5% 
increases) when considering only the MIPS 24 /zm 5cr 
detections. For visually identified disks and spheroids, 
the stellar mass increases by up to a factor of 2.7±0.4 
and 1.8±0.3, respectively. 

Ideally, one would like to disentangle what is the 



relative contribution to the size growth of a galaxy 
of newly-formed stars and system heating through 
merger/interactions. However, both processes arc prob- 
ably linked, since new star formation events are likely 
associated to the interactions that could inject energy to 
the systems. Consequently, a definitive answer to the 
problem of how galaxies grow require the help of elabo- 
rate modelling. It is worth saying, nonetheless, that both 
the observations at low and high-z show that galaxies 
have larger effective radii when observed in bluer bands. 
This means that younger stars arc preferentially located 
at larger galactocentric distances than older populations. 
In this paper, we have quantified how much the stellar 
mass grows through star formation events only. Once 
we reach a clear picture of how the galaxies can increase 
in size through mergers, our results will constrain the 
amount of stellar mass due to dry accretion that is neces- 
sary to migrate the high-z galaxies to the local size-mass 
relations. 

4. SUMMARY AND CONCLUSIONS 

We have analyzed the stellar mass growth in the form 
of newly-born stars in a sample of 831 if- band selected 
massive galaxies (A / (>10 11 M.q) as a function of struc- 
tural parameters (size and concentration). These galax- 
ies lie in the redshift range between z=0.2 and z^2. Our 
analysis is based on the measurement of the specific SFR 
for each galaxy based on their UV and IR emission, tak- 
ing advantage of the deep Spitzer data obtained by the 
FIDEL Spitzer/MlPS Legacy Project in the Extended 
Groth Strip. 

Our main results follow: 

— Most (more than 85% at any redshift) disk- 
like galaxies (identified by small Sersic indices, 
n<2.5) are detected by MIPS at 24 /im down to 
F(24)=15 ^Jy with a median flux F(24) = 190 fi3y. 

— A significant fraction (more than 55% at any 
redshift) of spheroid-like galaxies is detected at 
24 /im down to .F(24)=14 /iJy with a median flux 
F(24)=82 M Jy. 

— The MIPS detection fraction for spheroid-like 
galaxies is higher (70-90%) for larger (r e >5 kpc) 
galaxies, especially at z>l, where the detection 



fraction has a minimum around 30-40% for galax- 
ies with r e ^l kpc. No clear trend is found for disky 
galaxies of different sizes. 

There is basically no difference between the loci 
occupied by MIPS detected and undetected galax- 
ies in the stellar mass-size plane. However, for 
spheroid-like objects at a given stellar mass, MIPS 
non-detections are smaller than IR-bright sources 
by a factor of ^1.2. 

Most of the galaxies in our sample present spectral 
energy distributions which are consistent with an 
elliptical template in the UV/optical/NIR spectral 
range. Some galaxies morphologically classified as 
spheroids have UV emission tails which are typical 
of star- forming systems, most commonly at z>l. 

A ~10% fraction of the massive galaxies in our 
sample present X-ray or power-law-like mid-IR 
emission which must be linked to the presence of a 
bright (unobscured or obscured) AGN. 

Based on the measured specific SFRs, we estimate 
that spheroid-like galaxies have doubled (at the 
most, depending on the burst durations) their stel- 
lar mass due to newly-born stars alone from z~2 to 
z=0.2. Most of these mass increase (60%) occur at 
z>l, where specific SFRs are as high as 0.4 Gyr -1 . 

Disk-like galaxies have tripled (at the most) their 
stellar mass by newly- formed stars at z<2, with a 
more steady growth rate as a function of redshift. 
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